Magnetotunneling in a Two-Dimensional Electron-Hole System Near Equilibrium 
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We have measured the zero-bias differential tunneling conductance of InAs/AlSb/GaSb/AlSb/InAs 
heterostructures at low temperatures (1.7 K < T < 60 K) and under a magnetic field at various angles 
with the heterostructure's interfaces. Shubnikov-de Haas oscillations in the magnetoconductance 
reveal the two-dimensional (2D) character of the electrons accumulated at the InAs interfaces and 
yield their number in each of them. The temperature dependence of the oscillations suggests the 
formation of a field-induced energy gap at the Fermi level, similar to that observed before in simpler 
2D-2D tunneling systems. A calculation of the magnetoconductance that considers different 2D 
densities in the two InAs electrodes agrees with the main observations, but fails to explain features 
that might be related to the presence of 2D holes in the GaSb region. 
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Heterostructures based on the InAs-AlSb-GaSb system 
are at the heart of low-dimensionality interband tunnel- 
ing diodes. This system is also attractive from the point 
of view of basic physics because the unusual ordering of 
the InAs and GaSb energy bands allows to study the tun- 
neling properties of combined two-dimensional electron 
and hole gases (2DEG and 2DHG, respectively). 

Two alternative arrangements are usually consid- 
ered: GaSb/AlSb/InAs/AlSb/GaSb (type A) and 
InAs/AlSb/GaSb/AlSb/InAs (type B), in both of which 
electrons and holes accumulate at the InAs and GaSb in- 
terfaces, respectively. In type A, when the InAs layer is 
thin, a quasi-square InAs quantum well (QW) for elec- 
trons is created between AlSb barriers, and GaSb elec- 
trodes behave as sources of holes. In type B, a thin 
GaSb central layer forms a quasi-square QW for holes, 
whereas the InAs end layers act as electron sources. In ei- 
ther combination, when a voltage (V) is applied between 
the electrodes, negative differential resistance (NDR) ap- 
pears when the bandgap of GaSb starts blocking the tun- 
neling of electrcffls, from the conduction band of InAs into 
adjacent layers .§□ 

Until now, the emphasis has been on the study of 
the current-voltage (I-V) characteristics, often in con- 
nection with the materials properties of the hetesostruc- 
tures or with the details of their band structureJ3~@ The 
wealth of information available on the InAs-AlSb-GaSb 
system under an external bias constrasts with the limited 
knowledge about the details of the tunneling mechanism 
when the system is in quasi-equilibrium, that is, in the 
zero- voltage limit. Since carriers accumulate in quasi- 
triangular wells at the interfaces, the fact that tunnel- 
ing ideally occurs between quantized 2D states at the 
electmdes should be reflected in the conductance behav- 
ior .Q'E2l Moreover, since in both types of heterostructures 
vertical transport involves confined valence band states 
(holes in GaSb) and conduction band states (electrons 
in InAs), one could envision tunneling taking place via 
a 2DH(E)G -> 2DE(H)G -^ 2DH(E)G mechanism, with 
its own signature. 

The 2D nature of the electronic states at the central 



region, whether it is InAs or GaSb, is unquestioned, but 
the small confinement energy at the electrodes in type 
A heterostructures reduces (if not eliminates altogether) 
the 2D character of the holes U In type B structures, 
however, electrons in the quasi-triangular accumulation 
layers at the InAs electrodes are much more confined, and 
in a magnetic field their cyclotron energy is large enough 
to yield a distinct zero-dimensionality Landau-level spec- 
trum. These properties make InAs-GaSb-InAs based 
heterojunctions suitable to explore 2DEG -^ 2DIIG -^ 
2DEG tunneling, which is the objective of this work. 

We have measured the zero-bias differential magne- 
toconductance (G) of InAs/AlSb/GaSb/AlSb/InAs het- 
erostructures grown by molecular-beam epitaxy on p'^ 
GaSb substrates. The thickness of the undoped GaSb 
layer was 60 A and that of each AlSb barrier was 34 
A. Each electrode consisted of an 800 A thick undoped 
InAs region (residual n~ type doping estimated to be 1 
to 3 X lO^^cTTi"'^) adjacent to AlSb, followed by a thick 
(3000 A near the surface, 1 fim next to the substrate) n+ 
InAs (iVd = 2 X lO^^cm-3) layer. 

Figure 1 shows the potential energy profile of the het- 
erostructure along the growth direction, z. The overlap 
between the top of the GaSb valence band and the bot- 
tom of the InAs conduction band produces a self-limiting 
transfer of electrons from GaSb to InAs. As a result, 
an accumulation layer is formed at each InAs/ AlSb in- 
terface and holes are left behind in the GaSb quantum 
well. Ideally the total number of holes is twice that of 
electrons per interface. The number of transferred elec- 
trons increases with increasing overlap energy and with 
decreasing barrier thickness. In the figure, Eq indicates 
the ground-state energy for 2D electrons, Ep is the Fermi 
level of the system, and Hq represents the heavy-hole 
ground-state energy of the GaSb quantum well. The en- 
ergy values have been calculated by solving Poisson and 
Schrodinger equations self-consistently. 

The best evidence for the presence of holes in the cen- 
tral well is the observation (even at room -temperature) 
of NDR features in the I-V characteristics. □ As the volt- 
age between electrodes is increased, electrons from the 



emitter tunnel to the collector through the two barriers 
via the small energy window defined by Ep and if o ■ The 
tunneling process continues until the voltage is such that 
Ho is below the emitter's Eq level, at which point the 
tunnel current ceases and NDR appears. The shape of 
the I-V curves and the current-onset voltage depend on 
whether the in-plane momentum, fc|| , is conserved in the 
tunneling process. A detailed discussion of the similari- 
ties and differences between 2DEG -^ 2DEG and 2DEG 
-^ 2DHG is postponed for a future pubhcation. However, 
it is important to mention now that, in sharp contrast 
with the former case, tuncUing of momentum-conserving 
2D electrons into 2D hole states is possible for all volt- 
ages in the range between the threshold (which depends 
on the difference between the number of 2D electrons and 
holes) and NDR voltages. 
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FIG. 1. Energy-band profile of a type B heterostructure 
in the vicinity of the Fermi level Ef- The thickness of each 
AlSb barrier is 34 A. Eq and Hq represent the confined levels 
for two-dimensional (2D) electrons in the InAs conduction 
band (C.B.) and 2D holes in the GaSb valence band (V.B.), 
respectively. 



We focus here on the quasi-equilibrium conductance, 
that is, the slope of the TV characteristics in the V ^ 
limit. The measurements were done at low temperature 
(1.7K < T < 60K) by applying a small ac modulation 
voltage (200 [iV) between the electrodes of mesa-defined 
diodes with dimensions of 100 [im x 100 /iTTi, and a mag- 
netic field B, up to 11 T, forming an angle Lp with the 
plane of the interfaces. 

In Fig. 2(a) we show G(B) for (^ = 90° up to B = 11 T. 
The oscillations in G(B), discernible even below 0.5T, ex- 
hibit the characteristic B~^ periodicity of Shubnikov-de 
Haas oscillations. They resemble those previously ob- 
served in similar measurements on type A heterostruc- 
tures, which were explained in terms oLthc Landau lev- 
els formed in the InAs quantum well.O However, there 
are some notable differences between the oscillations ob- 



served in both types of heterostructures. While at low 
magnetic fields (B < 1 T) the amplitude of the oscilla- 
tions in type A increases monotonically with increasing 
field, that of the oscillations in Fig. 2(a) is modulated in 
a way similar to that observable in a two-subband sys- 
tem. Besides, a second set of oscillations is apparent in 
Fig. 2(a) at intermediate fields (IT < B < 4T), quite 
different from the spin-splitting doublets that appear at 
larger fields in type A. Finally, the oscillations in Fig. 
2(a) have an unusual shape at high B (B > 4T). 

The temperature and angular dependence of the low- 
field oscillations in Fig. 2(a) have elucidated their origin. 
The temperature dependence of the oscillations' ampli- 
tudes (see Fig. 3) yields effective mass values between 
0.025-0.028 m-o (?tio is the free-electron mass). The an- 
gular dependence reveals that the positions of the oscil- 
lations extrema are determined not by the total mag- 
netic field but by its perpendicular component i?j^, thus 
proving the 2D character of the charges involved in mag- 
nctotunneling. The low-field amplitude modulation and 
the doublets at intermediate fields indicate two differ- 
ent Landau-level ladders, each of them with its own B~^ 
periodicity and associated 2D carrier density. From the 
periods of the two oscillations, a simple estimate yields 
TVi = 5.13 X 10"cm-2 and A^z = 5.07 x lO^cm-^. 
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FIG. 2. Differential tunneling conductance in the 
zero-voltage limit (V=0) it vs a magnetic field perpendicular 
to the layers, for the heterostructure to which Fig.l. refers. 
(a) Experimental curve measured at T = 1.7 K. (b) Calculated 
curve using Eq. (1) in the text. The indices of the Landau 
levels in the two InAs electrodes are indicated as n and n'. 

These results are consistent with the magneto- 
conductance oscillations being due to tunneling between 
Landau levels of the 2D electron gases formed at the InAs 
electrodes. In this picture A^i and Ni represent the 2D 
electron densities of those electrodes; the small difference 
between the two indicates that the two interfaces are al- 



most identical and, therefore, that the subband energies, 
-Boi and i?02, are almost the same. Instead, it would be 
possible to interpret Ni and N2 as the carrier densities 
of two electron subbands of the same electrode, let us say 
the emitter. But if this were the case, the two densities 
should be quite different; indeed, our self-consistent cal- 
culations indicate that a second subband would be only 
marginally occupied. 

The above picture — 2D electrons tunneling from oc- 
cupied Landau levels in the emitter to equivalent empty 
states in the collector — is supported by the agreement 
found between the experimental results of Fig. 2(a) and 
the calculated tunneling magnetoconductance shown in. 
Fig. 2(b). The calculation is similar to that by LyoO 
for the conductance between two 2D electron gases in a 
perpendicular field B±, 



Gz2{B±) cc g{Bj_) V 6l„i^„2(")^P«i(M)Pn2(Ai- A), (1) 



where g{B±) — eB±/Trh is the Landau- level degeneracy; 
^ni.n2(^) is an "overlap" function between Landau levels 
ni and n2 of the electrodes, which, if the Landau-level 
index is conserved, is 1 when m = 712 and zero otherwise; 
p„. is the density of states (DOS) for the Landau level 
of index Ui ; and A is the difference between the subband 
energies of the two electrodes (A — \Eqi — Eq2\). We have 
used the approximation that only Landau levels at the 
Fermi energy fi contribute significantly to Gzz{B±), thus 
reducing the double sum on ni and n2 to Tii — nip and 
«2 — "2F . We have also assumed that the value of the 
overlap function is 1 because since the carrier densities 
in the two InAs layers are almost the same, in practice, 
at the Fermi energy the indices of the Landau levels of 
the two electrodes are the same. For the DOS we have 
used a Gaussian distribution, 
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where e„ — hujc{n + 1/2) is the energy of the n Landau 

level and F is the level damping (or broadening), given by 

1/2 
F = Q.hrjfiLOc/ B^ {rj is a sample-dependent parameter). 

The Fermi energy, which depends on magnetic field, was 

determined from the condition that the total number of 

carriers, 2A'o = A^i -f A^2, be constant: 
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where f{E,fj,,T) is the Fermi distribution function, and 
D^^{E) is the total density of states, which includes the 
DOS of each accumulation layer. Since the experimental 
results used for comparison were obtained at 1.7 K, in 
calculating fj, from Eq. (3) we assumed that T = OK. 

As seen in Fig. 2(b), this simple model reproduces well 
the peak positions throughout the entire field range. Our 
intention has not been to find systematically a "best" 



set of parameters but rather to illustrate how a set of 
reasonable values can explain the observed magnetocon- 
ductance oscillations. The following set has been used in 
the calculation of Fig. 2(b): 0.027 mo for the electron 
effective mass, 77 = 0.8 T^/^ for the parameter related 
with the broadening F, A^o = 5.7 x lO^^cm^^ and A = 
4 meV for the missalignement between the two energy 
subbands, Eqi and Eq2. The value used for the effective 
mass is consistent with that determined from the tem- 
perature dependence of the low-field Shubnikov-de Haas 
oscillations. The subband energy difference A is less than 
10% of any of the two energies; in other words, the differ- 
ence between A^i and N2 is quite small when compared 
to the average density in each electrode. The values of 
this difference and of the sum (2A''o), although somewhat 
larger than those estimated from the Shubnikov-de Haas 
periods, were chosen so that the calculated G(B) "fitted" 
well the experimental G(B) for the entire magnetic-field 
range. The small discrepancy may be due to our assump- 
tion that A'o is independent of B, an assumption that is 
especially questionable at high fields, given the physical 
origin of the 2D electrons (the transfer from the GaSb 
central well). 
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FIG. 3. Temperature dependence of the magnetoconduc- 
tance, relative to its zero-field value, for the heterostructure 
shown in Fig. 1. For clarity, the origin is different for each 
curve. While the amplitude of the low-field (below 3T) oscil- 
lations decreases monotonically with increasing T, that of the 
high-field (above 4T) broader oscillations first increases and 
then decreases. The features labeled as A, B and C, which 
are more distinguishable at the lowest T, cannot be explained 
by the simple model described in the text and may be related 
to the 2D holes in the GaSb central well. 



The different carrier densities of the two (nominally 



identical) accumulation layers suggest an asymmetry in 
the two AlSb barriers, a fact confirmed by the asymme- 
try observed between the two polarities of the diodes' 
I-V characteristics. If one barrier were thicker or higher 
than the other, then the intrinsic transfer of electrons 
from GaSb to each InAs interface would be different. In 
practice, barrier asymmetry can happen for a number of 
reasons; for instance, Ga segregation during the epitaxial 
growth of the heterostructures would lead to an Ga-rich 
AlGaSb layer, which has a smaller bandgap than AlSb 
and would present a smaller tunneling barrier. This seg- 
regation process is known to be the more significant the 
higher the substrate temperatureB Since our heterostruc- 
tures were prepared at 500 °C — a relatively high temper- 
ature for these materials — it would not be surprising that 
the barrier asymmetry is caused by that mechanism. 

In spite of the overall good agreement between the ex- 
periment and the calculation, at high magnetic fields the 
measured conductance behaves in a way for which the 
model cannot account. In Fig. 2(a) there are three peaks 
(marked by "A" , "B" , and "C" ) that do not have a coun- 
terpart in Fig. 2(b). Although the fact that they were 
clearly resolved only at the lowest temperatures (see Fig. 
3) might suggest spin splitting effects, the inclusion of 
spin in the calculation (not shown in Fig. 2(b)) failed to 
explain both the number of peaks and the fields at which 
they occur. 

The temperature dependence of the broader Shubnikov 
-de Haas oscillations underlying peaks A, B, and C also 
warrants a comment. As T increases, the amplitude of 
the oscillations increases a bit before it decreases at high 
temperatures (see Fig. 3), suggesting a thermally acti- 
vated behavior followed by the usual T-induced weaken- 
ing of quantum effects, when fcsT becomes comparable to 
the cyclotron energy and the scattering time is compara- 
ble to the cyclotron time. Although the low-temperature 
range investigated here is too limited to draw any conclu- 
sion, it is quite possible that the activated behavior is a 
consequence of the fomation of a field-induced energy gap 
at the Fermi level. Such a gap has been proposed to ex- 
plain a similar temperature dependence of the 2DEG p^ 
2DEG conductance in GaAs-GaAlAs heterostructures .EJ 
Calculations have shown that the gap appears as a result 
of the additional energy cost required to extract an elec- 
tron from the emitter Fermi sea (creating a "vacancy" ) 
and to inject it into the collector one (creating an "inter- 
stitial") when many-body effects (for example, Coulomb 
interactions) become sizeable.Ej 

The quasi-equilibrium magnetoconductance experi- 
ments presented here show unambiguously the partici- 
pation of 2D electrons in the 2DEG -^ 2DHG -^ 2DEG 
tunneling process, but do not reveal the role of the 2D 
holes in it. Whether the unaccounted low-temperature 
peaks (A, B and C) are evidence of the holes we do not 
know. However, it is worth pointing that, based on their 
shape and the temperature range in which they are most 
visible, those peaks resemble the features found in the in- 
plane magnetoresistance of InAs-GaSb heterostructures 



in which 2D electrons and holes coexist .113 Although these 
features remain unexplained, their connection with 2D 
holes is well established, since they disappear when the 
2D holes [are eliminated by the application of hydrostatic 
pressure.tJ It is then reasonable to imagine that the ex- 
tra peaks in Fig. 3 are somehow related to the quantized 
holes in the GaSb quantum well, but more work, both 
theoretical and experimental, is needed to set this possi- 
ble correspondence on a firmer ground. 



ACKNOWLEDGMENTS 

We thank Dr. A. Sacedon for her participation in the 
early stages of this work, Dr.W. I. Wang for contributing 
the heterostructures used in it, and Dr. Aleiner for useful 
discussions. The support of the US Army Research Office 
and of a NATO postdoctoral fellowship (to E. M. G.) 
have made this work possible. 



Permanent address: Universidad Complutense, Madrid, 
Spain. 
^ J. R. Soderstrom, D. H. Chow, and T. C. McGill, Appl. 

Phys. Lett. 55, 1094 (1989). 
^ L. F. Luo, R. Beresford, and W. I. Wang, Appl. Phys. Lett. 

55, 2023 (1989). 
^ B. Z. Nosho ei al, J. Vac. Sci. Technol. B 16, 2381 (1998). 
* H. Kitabayashi, T. Waho, and M. Yamamoto, J. Appl. 

Phys. 84, 1460 (1998). 
5 M. Ogawa, T. Sugano, and T. Miyoshi, Physica E 7, 840 

(2000). 
''T. Takamasu et al., Physica B 184, 259 (1993); 201, 380 

(1993). 
^ R. R. Marquardt et al, Phys. Rev. B 53, 13624 (1996); Y. 

X. Liu et al, Phys. Rev. B 55, 7073 (1997). 
® E. E. Mendez et al, Physica E 6, 335 (2000). 
^ E. E. Mendez and L. L. Chang, Surf. Sci. 229, 173 (1990). 
^° J. P. Eisenstein et al, Phys. Rev. B 44, 6511 (1991). 
" E. E. Mendez, J. Nocera, and W. I. Wang, Phys. Rev. B 

45, 3910 (1992). 
^^E. E. Mendez et al, Phys. Rev. B 43, 5196 (1991); E. E. 
Mendez, J. Nocera, and W. I. Wang, Phys. Rev. B 47, 
13937 (1993). 
" S. K. Lyo, Phys. Rev. B 57, 9114 (1998). 
" J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, Phys. Rev. 

Lett. 69, 3804 (1992). 
^^ S.-R. E. Yang and A. H. MacDonald, Phys. Rev. Lett. 70, 

4110 (1993). 
^® E. E. Mendez, L. Esaki, and L. L. Chang, Phys. Rev. Lett. 

55, 2216 (1985). 
^■^ J. Beerens et al, Phys. Rev. B 35, 3039 (1987). 



